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It is shown that quark mass dependence induced by one loop corrections to the Breit-Fermi spin- 
dependent one gluon exchange potential permits an accurate determination of heavy-light meson 
masses. Thus the D s (2317) is a canonical cs meson in this scenario. The multiplet splitting rela- 
tionship of chiral doublet models, M(l + ) — M(l~) = M(0 + ) — M(0~), holds to good accuracy in 
the D and D 3 systems, but is accidental. Radiative transitions and bottom flavoured meson masses 
are discussed. 



I. INTRODUCTION 



o ■ 
o : 

BaBar's discovery of the D s (2317) statepj generated strong interest in heavy meson spectroscopy chiefly due to its 
surprisingly low mass with respect to expectations. These expectations are based on quark models or lattice gauge 
theory. Unfortunately, at present large lattice systematic errors do not allow a determination of the D s mass with 
a precision better than several hundred MeV. And, although quark models appear to be exceptionally accurate in 
describing charmonia, they are less constrained by experiment and on a weaker theoretical footing in the open charm 
sector. It is therefore imperative to examine reasonable alternative descriptions of the open charm sector. 

The -D s (2317) was produced in e + e _ scattering and discovered in the isospin violating D s tt decay mode in KKttt: 
and KKititit mass distributions. Its width is less than 10 MeV and it is likely that the quantum numbers are 
J p = O + 0- Finally, if the D s tt° mode dominates the width of the D s (2317) then the measured product of branching 



' ratios 01 

00 : 

Br(B° -> D S (2317)K) ■ Br{D s (2317) -> D s n°) = (4.4 ± 0.8 ± 1.1) • 10~ 5 (1) 

implies that Br(B — > 15,5(2317)^) w Br(B — > D S K), consistent with the D s (2317) being a canonical + cs meson. 
In view of this, Cahn and Jackson have examined the feasibility of describing the masses and decay widths of the 
Oh, low lying D and D s states within the constituent quark modelQ. They assume a standard spin-dependent structure 
q^. for the quark-antiquark interaction (see below) and allow general vector and scalar potentials. Their conclusion is 
| that it is very difficult to describe the data in this scenario. 

• Indeed, the D s (2317) lies some 160 MeV below most model predictions (see Ref.[2j for a summary), leading to 
speculation that the state could be a DK molecule^ or a tetraquark[(J. Such speculation is supported by the isospin 
• rH , violating discovery mode of the -D s (2317) and the proximity of the S-wave DK threshold at 2358-2367 MeV. Other 
studies have been made with QCD sum rules0, using heavy quark symmetry to examine decay models^, or in 
unitarised chiral models 0. 

Although these proposals have several attractive features, it is important to exhaust possible canonical cs descrip- 
tions of the D s (2317) before resorting to more exotic models. Here we propose a simple modification to the standard 
vector Coulomb+scalar linear quark potential model that maintains good agreement with the charmonium spectrum 
and agrees remarkably well with the D and D s spectra. Possible experimental tests of this scenario are discussed. 



II. A QUARK MODEL OF OPEN CHARM STATES 

The spectra we seek to explain are summarised in Table [|] Unfortunately, the masses of the Do (labelled a) and D\ 
(labelled b) are poorly determined. Belle have observed 10J the Do in B decays and claim a mass of 2308 ± 17±32 MeV 
with a width of T = 276±21±18±60, while FOCUS[I3 nnd 2407±21±35 MeV with a width T = 240±55±59. While 
some authors choose to average these values, we regard them as incompatible and consider the cases separately below. 
Finally, there is an older mass determination from Belle fll of ' 2290±22±20 MeV with a width of T = 305±30±25. The 
D[ has been seen in B decays to Dtttt and D*tttt by Belle|l3j| . A Breit-Wigner fit yields a mass of 2427 ± 26 ± 20 ± 15 
MeV and a width of 384+qq 7 ± 24 ± 70 MeV. Alternatively, a preliminary report from CLEO[lJ| cites a mass of 
2461±34 ± 10 ± 32 MeV and a width of 290±2.^ ± 26 ± 36 MeV. Finally, FOCUS 15] obtain a lower neutral D[ mass 
of 2407 ± 21 ± 35 MeV. Other masses in Table |T| are obtained from the PDG compilation 0. 

In addition to the unexpectedly low mass of the D s (2317), the _D S (2460) is also somewhat below predictions 
assuming it is the D s i (Godfrey and Isgur, for example, predict a mass of 2530 MeV^T})- It is possible that an 
analogous situation holds in the D spectrum, depending on the mass of the Do- The quark model explanation of these 
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TABLE I: Low Lying D and D s Spectra (MeV) 

jPH 0~ I" 0+ 1+ 1+ 2+ 

D 1869.3 ±0.5 2010.0 ±0.5 a b 2422.2 ± 1.8 2459 ± 4 

D s 1968.5 ± 0.6 2112.4 ± 0.7 2317.4 ± 0.9 2459.3 ± 1.3 2535.35 ± 0.34 2572.4 ± 1.5 



states rests on P-wave mass splittings induced by spin-dependent interactions. A common model of spin-dependence 
is based on the Breit-Fermi reduction of the one-gluon-exchange interaction supplemented with the spin-dependence 
due to a scalar current confinement interaction. The general form of this potential has been computed by Eichten and 
Feinberg|l8j| at tree level using Wilson loop methodology. The result is parameterised in terms of four nonperturbative 
matrix elements, Vi, which can be determined by electric and magnetic field insertions on quark lines in the Wilson 
loop. Subsequently, Gupta and Radford 19] performed a one-loop computation of the heavy quark interaction and 
showed that a fifth interaction, V5 is present in the case of unequal quark masses. The net result is a quark-antiquark 
interaction that can be written as: 



Vqq = V conf ± V SD (2) 

where V con f is the standard Coulomb±linear scalar form: 

VW(r) = ~— + &r (3) 
o v 

and 
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Here L = L g = — Lq, r = |r 



Y q — Yq\ is the QQ separation and the Vi — Vi(m q , m^; r) are the Wilson loop matrix 
elements discussed above. 

The first four Vi are order a s in perturbation theory, while V5 is order a 2 s ; for this reason V5 has been largely 
ignored by quark modellers. The exceptions known to us are Ref. |20| . which examines S-wave masses for a variety 
of heavy-light mesons in a model very similar to that presented here, and the second of Ref. |19| . which does not 
consider scalar confinement contributions to the spin-dependent interaction. More recently, Cahn and Jacksonj^ only 
consider V\ - V4 in an analysis of the D s system. In practice this is acceptable (as we show below) except in the case 
of unequal quark masses, where the additional spin-orbit interaction can play an important role. 

Here we propose to take the spin-dependence of Eqn. 0] seriously and examine its effect on low-lying heavy-light 
mesons. Our model can be described in terms of vector and scalar kernels defined by 



V conf = V + S (5) 

where V = — 4a s /3r is the vector kernel and S = br is the scalar kernel, and by the order a 2 s contributions to the Vi, 
denoted by SVi. Expressions for the matrix elements of the spin-dependent interaction are then 



Vi = -S + SVi 

v 2 = v + sv 2 

V 3 = V'/r-V" + 5V 3 

Vi = 2V 2 V + 5V4 

V 5 = SV 5 . (6) 



Explicitly, 
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where Cf = 4/3, Ca = 3, &o = 9, 7 F = 0.5772, and the scale p, has been set to 1 GeV. 

The hyperfine interaction (proportional to V4) contains a delta function in configuration space and is normally 
'smeared' to make it nonperturbatively tractable. For this reason we choose not to include SV4 in the model definition 
of Eqn. In the following, the hyperfine interaction has been treated nonperturbatively and the remaining spin- 
dependent terms are evaluated in leading-order perturbation theory. 

We have confirmed that the additional features do not ruin previous agreement with the charmonium spectrum. 
For example, Ref. obtains very good agreement with experiment for parameters m c — 1.4794 GeV, a s = 0.5461, 
b = 0.1425 GeV 2 , and a = 1.0946 GeV. Employing the model of Eqn. [7| worsens the agreement with experiment, 
but the original good fit is recovered upon slightly modifying parameters (the refit parameters are m c = 1.57 GeV, 
a s = 0.52, b = 0.15 GeV 2 , and a = 1.3 GeV). 



TABLE II: Model Parameters 



model 


a s 


b (GeV) 


a (GeV) 


m c (GeV) C (GeV) 


low 


0.46 


0.145 


1.20 


1.40 -0.298 


avg 


0.50 


0.140 


1.17 


1.43 -0.275 


high 


0.53 


0.135 


1.13 


1.45 -0.254 



The low lying cs and cu states are fit reasonably well with the parameters labelled 'avg' in Table [H] Predicted 
masses are given in Table ITTT1 Parameters labelled 'low' in Table [D] fit the D mesons very well, whereas those labelled 
'high' fit the known D s mesons well. It is thus reassuring that these parameter sets are reasonably similar to each 
other and to the refit charmonium parameters. (Note that constant shifts in each flavour sector are determined by 
the relevant pseudoscalar masses.) 

The predicted D s q mass is 2341 MeV, 140 MeV lower than the prediction of Godfrey and Isgur and only 24 MeV 
higher than experiment. We remark that the best fit to the D spectrum predicts a mass of 2287 MeV for the D 
meson, in good agreement with the preliminary Belle measurement of 2290 MeV, 21 MeV below the current Belle 
mass, and in disagreement with the FOCUS mass of 2407 MeV. 

The average error in the predicted P-wave masses is less than 1%. It thus appears likely that a simple modification 
to the spin-dependent quark interaction is capable of describing heavy-light mesons with reasonable accuracy. 



TABLE III: Predicted Low Lying Charm Meson Spectra (GeV) 



flavour 


0~ 


1" 


0+ 


1+ 


1+ 


2+ 


D 


1.869 


2.017 


2.260 


2.406 


2.445 


2.493 


D a 


1.968 


2.105 


2.341 


2.475 


2.514 


2.563 



We examine the new model in more detail by computing P-wave meson masses (with respect to the ground state 
vector) as a function of the heavy quark mass. Results for Qu and Qs systems are displayed in Fig. ^ All panels 
indicate that the approach to the heavy quark limit is very slow. In the case of the traditional Qu system the heavy 
quark doublets are inverted (with the j q = 1/2 doublet higher than the j q — 3/2), in disagreement with experiment. 
Alternatively, the one-loop model displays the expected heavy quark behaviour. Furthermore, the predicted mass 
splittings at the charm quark scale are near experiment for D masses (points on the panels). A similar situation holds 
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for the D s system (right panels), except in this case it is the D s and D' s that do not agree with the traditional model 
predictions. 

Although the reliability of the model is suspect in the case of light Q masses, it is intriguing that the one-loop 
model scalar- vector mass difference gets small in this limit. Thus it is possible that the enigmatic ao and fo mesons 
may simply be qq states. 

Finally, one obtains M(h c ) > M(xci) in one-loop and traditional models, in agreement with experiment. However, 
experimentally M(/i) - M(hi) « 100 MeV and M(ai) - M(bx) w MeV, indicating that the 3 Pi state is heavier 
than (or nearly degenerate with) the 1 P\ light meson state. Thus the sign of the combination of tensor and spin-orbit 
terms that drives this splitting must change when going from charm quark to light quark masses. This change is 
approximately correctly reproduced in the traditional model (lower left panel of Fig. The one-loop model does 
not reproduce the desired cross over, although it does come close, and manipulating model parameters can probably 
reproduce this behaviour. We do not pursue this here since the focus is on heavy-light mesons. 




456789 10 012345 
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FIG. 1: M(P-wave) - M(vector) as a Function of the Heavy Quark Mass. D System (left); D s System (right); one-loop model 
(top); traditional model (bottom). D2 (dotted line); D[ (thin dashed line); Di (dashed line); Do (solid line). The points are 
experimental mass splittings. 



III. MIXING ANGLES AND RADIATIVE DECAYS 

The lack of charge conjugation symmetry implies that two nearby low lying axial vector states exist (generically 
denoted as D\ and D[ in the following). The mixing angle between these states can be computed and compared to 
experiment (with the help of additional model assumptions). We dehne the mixing angle via the relations: 

|Di) = +cos($| 1 P 1 > +sim»| 3 P 1 ) 

\D[) =-sin(0)| 1 F 1 )+cos( ( /))| 3 P 1 ). (8) 

In the following, we choose to define the D[ as the heavier axial state in the heavy quark limit. In this limit a 
particular mixing angle follows from the quark mass dependence of the spin-orbit and tensor terms, <Phq = —54.7° 
(35.3°), if the expectation of the heavy-quark spin-orbit interaction is positive (negative). It is often assumed that 
the heavy quark mixing angle holds for charmed mesons. 

Fig. |2]shows the dependence of the mixing angle on the heavy quark mass for Qu and Qs mesons for the traditional 
and extended models. The effect of the one-loop terms is dramatic: for the Qu system the relevant spin-orbit 
matrix element changes sign, causing the heavy quark limit to switch from 35.3° to —54.7°. Alternatively, both 
models approach —54.7° in the Qs system. There is strong deviation from the heavy quark limit in both cases: 
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4>(D S ) w 0(1?) ~ —70°. This result is not close to the heavy quark limit (which is approached very slowly) - indeed 
it is reasonably close to the unmixed limit of ±90°! 
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FIG. 2: Di (left) and D s \ (right) Mixing Angles vs. Heavy Quark Mass. The traditional model is given by the dashed line; 
the extended model is the solid line. 

Mixing angles can be measured with the aid of strong or radiative decays. For example, the D[ is a relatively 
narrow state, T(I? 1 ) = 20.4 ± 1.7 MeV, while the Di is v ery broad. This phenomenon is expected in the heavy quark 
limit of the 3 Pq and Cornell strong decay models 0, H3- Unfortunately, it is difficult to exploit these widths to 
measure the mixing angle because strong decay models are rather imprecise. 

Radiative decays are possibly more accurate probes of mixing angles because the decay vertex is established and the 
impulse approximation has a long history of success. Table llVl presents the results of two computations of radiative 
decays of D and D s mesons. Meson wavefunctions are computed with 'average' parameters, as above. Transition 
matrix elements arc evaluated in the impulse approximation and full recoil is allowed. The column labelled 'nonrel' 
reports transition matrix elements computed in the nonrelativistic limit, while the column labelled 'rel' contains 
results obtained with the full spinor structure at the photon vertex. 

The nonrelativistic results can differ substantially from those of Refs. jU because those computations were 
made in the zero recoil limit where an El transition, for example, is diagonal in spin. Thus the decay D\ —* D*^ can 
only proceed via the 3 Pi component of the D\. Alternatively, the computations made here are at nonzero recoil and 
hence permit both components of the D\ to contribute to this decay. The table entries indicate that nonzero recoil 
effects can be surprisingly large. 



TABLE IV: Meson Radiative Decay rates (keV). 



mode 


g 7 (MeV) 


nonrel 


rel 


expt 


D + * -> 7D+ 


136 


1.38 


0.08 


1.5 ±0.5 


D°* -» jD° 


137 


32.2 


13.3 


< 800 


D+ 1 D*+ 


361 


76.0 


7.55 




D° -» 1 D*° 


326 


1182 


506. 




D+ -» 7 £>*+ 


381 


(6.34s) 2 + (3.22s + 5.9c) 2 
(27.05s) 2 + (19.33s + 9.63c) 2 
(6.34c) 2 + (-3.22c + 5.9s) 2 
(27.26c) 2 + (19.35c -9.83s) 2 


(2.00s - 0.13c) 2 + (0.13s + 4.23c) 2 
(17.65s - 0.15c) 2 + (12.28s + 6.01c) 2 
(2.00c + 0.13s) 2 + (-0.13c + 4.23s) 2 
(17.78c + 0.15s) 2 + (12.29c - 6.13s) 2 




D° -* 1 D*° 


380 




D '+ _ 7 zr+ 


381 




D '0 _ lD *o 


384 




D+ -> jD+ 


494 


(5.49s + 4.75c) 2 


(4.17s -0.60c) 2 
(5.56s + 18.78c) 2 
(4.17c -0.60s) 2 




D° -» 1 D° 


493 


(8.78s + 31.42c) 2 




D'+ -* 7 L»+ 


494 


(-5.49c + 4.75s) 2 




Df 7 D° 


498 


(-8.90c + 31.41s) 2 


(-5.62c + 18.78s) 2 




D+ -> 7 D*+ 


413 


15.0 


6.49 




D° 2 -» 1 D*° 


412 


517 


206 




D* s -» jD a 


139 


0.20 


0.00 




Dso -> jD* s 


196 


6.85 


0.16 




Dsi -> ~/D* s 


322 


(1.84s) 2 + (0.99s + 2.39c) 2 
(2.13c) 2 + (-0.87c + 3.62s) 2 


(0.18s - 0.07c) 2 + (-0.44s + 2.13c) 2 
(0.24c - 0.10s) 2 + (0.64c + 3.19s) 2 




D' sl jD2 


388 




Dsi — ► 7-Ds 


441 


(2.68s + 1.37c) 2 
(3.54c - 1.12s) 2 


(2.55s - 1.21c) 2 




D' sl ^ -/D s 


503 


(3.33c + 1.52s) 2 




Ds2 iD* s 


420 


1.98 


3.94 
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Further complicating the analysis is the large difference seen between the nonrelativistic and relativistic models 
(see, eg, D + * — > ^D + ). This unfortunate circumstance is due to differing signs between the heavy and light quark 
impulse approximation subamplitudes. Employing the full quark spinors leaves the heavy quark subamplitude largely 
unchanged, whereas the light quark subamplitude becomes larger, thereby reducing the full amplitude. The effect 
appears to be at odds with the only available experimental datum (D* — ► Dj) . Clearly it would be useful to measure 
as many radiative transitions as possible in these sectors to better evaluate the efficacy of these (and other) models. 

Once the decay model reliability has been established, ratios such as T(Di — > jD*)/r(D[ — > jD*) and T(Di — > 
7-D) /T(D\ — > 7Z?) will help determine the D± mixing angle. 



IV. DISCUSSION AND CONCLUSIONS 



A popular model of the D s mesons is based on an effective lagrangian description of mesonic fields in the chiral and 
heavy quark limits|25j. Deviations from these limits induce mass splittings which imply that the axial-vector and 
scalar-pseudoscalar mass differences are the same. Since the premise of this idea has been questioned in Refs. 
it is of interest to consider this mass difference in the present model. Splittings for the three parameter sets considered 
above are shown in Table [V] Evidently, the chiral multiplet relationship holds to a very good approximation in both 
the D and D s sectors and is robust against variations in the model parameters. 

Nevertheless, the near equivalence of these mass differences must be regarded as an accident. Indeed, the B masses 
given in Table IVTl indicate that this relationship no longer holds. It would thus be of interest to find P-wave open 
bottom mesons (especially scalars). These data will distinguish chiral multiplet models from the model presented 
here and from more traditional constituent quark models. For example, Godfrey and Isgur claim that the Bo meson 
lies between 5760 and 5800 MeV; the B sQ mass is 5840-5880 MeV, and the B c0 mass is 6730-6770 MeV. Of these, our 
B s q mass is predicted to be 65-105 MeV lower than the Godfrey-Isgur mass. 



TABLE V: Chiral Multiplet Splittings (MeV). 



params 


M(l+(l/2+))-M(l-; 


M(0+) -M(0~) 


D low 


411 


412 


D avg 


391 


389 


D high 


366 


368 


D a low 


384 


380 


D a avg 


373 


370 


D a high 


349 


346 



TABLE VI: Low Lying Bottom Meson Masses (MeV) 



flavour 


0" 


l - 


0+ 


1+ 


1+ 


2+ 


B 


5279 


5322 


5730 


5752 


5753 


5759 


expt 


5279 


5325 




5724 ± 4 ± 7 




5748 ± 12 


B a 


5370 


5416 


5776 


5803 


5843 


5852 


expt 


5369.6 


5416.6 










B c 


6286 


6333 


6711 


6746 


6781 


6797 


expt 


6286 













The bottom flavoured meson spectra of Table Ivfl have been obtained with the 'average' extended model parameters 
and TOfc = 4.98 GeV. As with the open charm spectra, a flavour-dependent constant was fit to each pseudoscalar. 
The second row reports recently measured P-wave B meson masses 27]; these are in reasonable agreement with the 
predictions of the first row. 

When these results are (perhaps incorrectly) extrapolated to light quark masses, light scalar mesons are possible. 
Thus a simple qq interpretation of the enigmatic ao and fo mesons becomes feasible. 

Finally, the work presented here may explain the difficulty in accurately computing the mass of the D s o in lattice 
simulations. If the extended quark model is correct, it implies that important mass and spin-dependent interactions 
are present in the one-loop level one-gluon-exchange quark interaction. It is possible that current lattice computations 
are not sufficiently sensitive to the ultraviolet behaviour of QCD to capture this physics. The problem is exacerbated 
by the nearby, and presumably strongly coupled, DK continuum; which requires simulations sensitive to the infrared 
behaviour of QCD. Thus heavy-light mesons probe a range of QCD scales and make an ideal laboratory for improving 
our understanding of the strong interaction. 
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